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ABSTRACT 

We explore thermal X-ray iron line emission from the galactic X-ray binary GX 339-4 in 
the off state, using the models of the advection dominated accretion flow (ADAF) without or 
with outflows. The equivalent widths of hydrogen-like and helium-like thermal iron lines are 
calculated with different model parameters including viscosity parameter a, mass accretion rate 
at the outer radius of the ADAF mout and outflow strength parameter p. Our calculations show 
that the equivalent widths of thermal iron lines emitted from the pure ADAF, i.e., the ADAF 
without outflows, should be very small, assuming a solar metallicity for the accreting gas in 
the accretion flow. Strong thermal iron lines are expected to be emitted from the ADAF with 
relatively strong outflows. For reasonable choice of parameters, the total equivalent width of the 
He-like and H-likc thermal iron lines reaches to > 500 eV for accreting gas with solar metallicity. 
The observation of strong thermal X-ray lines from GX 339-4 at the off state may give a clue 
to the accretion mode of the source and provide evidence for the presence of outflows/winds in 
the accretion flow around the black hole in GX 339-4. It is found that the values of rhout and p 
are degenerate, i.e., the observed X-ray continuum spectrum can be fairly well reproduced with 
different sets of the parameters rhout and p. Such degeneracy can be broken when the thermal 
X-ray line emission data is available. We also compare our results with those in the previous 
similar work. 
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The galactic X-ray binary (XRB) GX 339- 
4, one of the earliest proposed black hole can- 
didates, has been attracted much attention af- 
ter discovered in 19 71 by the OSO-7 satellite 
([Markert et al.l Il973f) . Based on the long-term 
broadband energy spectral observations, from 
the radio, infrared and optical, to the X-rays, 



ous issues ( 


Fender et al. 1999: 


Corbel et al.l 20001: 


Kong et al. 


2000l iNowak et al. 


20021: iMiller et al. 


2004. 2006: 


Dunn et al.l2008; Tomsick et al.ll2008). 



state, but it spends most of its time in the 
low/hard state. A common argument is that 
the spectral states transitions are related to the 
changes of the mass accretion rate and the ge- 
ometry of the accretion flow, though the physi- 
cal mech anism of the changes are still quite un- 
certain dNaravan. Mahadevan fc Quataerd Il998t 
Donell200l[ ). 

Many different models were suggested to ex- 
plain the obs ervations of the different spec- 
tral states ( see iNaravan. Mahadevan fc Quataert 



Its spectral and temporal X-ray properties show 
a wide variety of canonical XRB states, transit- 
ing from off/quiescent state, low/hard state, in- 
termediate state, high/soft state, and very high 



1998 : Done! 2001 . and references therein). The 
disk-|-corona models with different geometries, 
such as, an accretion disk corona (ADC) above 
or surrounded by a cold thin disk , or ot her sim- 
ilar geometries (see Nowak et al.l l2002l for de- 
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tails and references therein), are used to interpret 
the observational features of the accretion flows 
around the black hole in different states. The ob- 
served X-ray spectra in different spectral states 
can be explained by the Comptonization ( e.g., 
soft photons from the cold thin accretion disk are 
Comptonized in a hot corona) wit h different model 



parameters (INowak et al.l 120021 : ICorongiu et al 



2003h . iNaravanI (Il996[ r proposed that these differ 



ent spectral states can be explained by outer thin 
disk plus inner advection dominated accretion flow 
(ADAF) models with varying acc r etion rate and 
transition radius (jEsin et al.l ll997L 119981 ) . In this 
scenario, the transition radius increases with de- 
creasing mass accretion rate m {ifi = M/M-Edd)- 
There is a critical accretion rate merit, above 
which the ADAF is suppressed and a standard 
thin disk is p resent in the inner region of the ac- 
cretion flow. iMiller et al.l (|2004l . |2006[ ) suggested 
that a standard thin disk may extend at or near 
to the innermost stable circular orbit (ISCO), at 
least in bright phases of the low/hard state, af- 
ter fitting the spectra of GX 339-4 observed by 
Chandra during its decline of 2002-2003 outburst 
and by XMM-Newton during its 2004 outburst 
with a cool thin disk component and a relativis- 
tic Fe K line model. It is shown that the disk 
reflection model can also fit the spectra in the 
hard state when the total 1-100 keV luni inosity 
even lower than 0.023LEdd (|Tomsick et al JHoos). 
Moreover, the presence of jet/outfiow in the accre- 
tion fiows has been explored to explain the radio 
emission and th e radio/X-ray corre l ations in the 
low/hard state ( Fender et al. _ 1999 Corbel et al 



20001 : iMarkoff et al.ll2003l : ICorbel et al.ll2003h 

Fluorescent Fe Ka line emission is a very pow- 
erful trace of the inner edge (i?in) of the thin accre- 
tion disk around the black holes. iTomsick et al 
(,200 9) used Suzaku and RXTE observations of 
GX 339-4 during 2008 to study the iron line pro- 
file at low luminosity state when the total 1-100 
keV luminosity is 0.0014LEdd, and showed that 
the inner edge of the thin disk recedes and the 
truncated accretion disk is present around a black 
hole at low luminosity. If the ADAF exists in 
th e inner region of the accretion flow as suggested 
bv iNaravan (1996), thermal X-ray iron lines may 
also emitted from the hot plasma in the ADAF 
surroun ding the black hole at very low lumi- 
nosity ( Naravan fc Raymond 19991 Perna et al 



20001: IXu et al.l l2006t IXu fc Cad l2009l) . It is 

shown that the average 0.4-12 keV X-ray flux 
that measured by RXTE between 2001 March 
and 2002 February and also on 2003 September 
29 is - 5.4 X lO-i^erg cm'^ s-\ when GX 339-4 



was in the quiescent state ( Tomsick et al. 2009( ). 
This average flux corresponds to a luminosity of 
Lx ^ 4.2 X lO^^erg s^^ - 3.2 x lO'^Ledd, if we 
adopt values of m = 10 (m = M/Mq) and d — S 
kpc for GX 339-4. Chandra even detected a 0.4- 
11 keV flux of 3.8 x lO'^erg cm -^ s'^ on 2003 



September 29 (jGallo et al.l 120031 ). which is more 
than an order of magnitude lower than the average 
quiescent value above. In such very low luminosity 
quiescent states, the accretion rate is much lower 
than the critical accretion rate TTicrit, the outer 
thin disk may be suppressed, and the temperature 
of the plasma in the outer region of the ADAF 
may be low to ~ lO^'^K, so the He-like (Fe 
XXV) and H-like (Fe XXVI) thermal X-ray iron 
lines centered at energy ~ 6.7 keV and ~ 6.97 keV, 
respecti vely, can be produced i n this r egion of the 



ADAF. iNaravan fc Ravmondl (|l999h calculated 



the thermal X-ray line emission including He-like 
and H-like iron lines from another X-ray binary, 
V404 Cyg, in quiesecence with the ADAF models 
with and without winds, and concluded that an 
ADAF with a large outer radius would produce 
reasonably strong thermal iron lines if winds are 
present in the source. However, for the ADAF 
models with winds, the dynamical structure of 
the ADAFs ar e solved based on some simp lified 
assumptions in Naravan fc Ravmond I ( 19991) . In 
this work, we will use the improved numerical 
approach to derive the global two-temperature 
structures of the ADAF models by solving the 
full relativistic hydrodynamical equat ions with 
themodeling techniques suggested by iManmotol 
(|2000f ). in which all the radiation processes are 
included. 

In this work, we model the thermal X-ray iron 
line emission from the galactic X-ray binary GX 
339-4 in the off state, using the advection domi- 
nated accreti on flow (ADAF) m odels with or with- 
out outflows (jBlandford fc Beg clman 1999). The 
models are constrained by the Chandra's observed 
X-ray continuum, the power-law photon index of 
2 ± 0.2 and the integrat ed 0.4-11 keV fl ux of 



3.8 x 10-"erg cm-^ g-i (IGallo et al.l 120031) . We 



describe the calculation of dynamical structures of 
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the ADAF in §2. The calculation of the equivalent 
widths of the thermal X-ray iron line emission is 
described in §3. The results are shown in §4. In 
§5., we discuss the results and the physical im- 
plications of the results, a nd compare our results 



of V4 04 Cyg with those in iNaravan fc Raymond 
(|l999h . 



2. Dynamical structures of the ADAF 
models 



We employ the approach suggested by Manmotol 
(|2nno h to calculate the global structure of an ac- 
cretion flow surrounding a Schwarzschild black 
hole in the general relativistic frame. All the radi- 
ation processes are included in the calc ulations of 
the g l obal accretion flo w structure (see lManmot"ol 



tne g lob 
20001: 1^ 



"uan et al. I l2009l for details and the refer- 
In this work, we include outflows 

which 



ences therein) 

in the calculation of the ADA F structure 

have not been considered by lManmot"ol (|2000t ). 
and the values of some paramet ers adopted in this 
work are different from those in iManmotol ([2000, ). 

The global structure of an accretion flow sur- 
rounding a black hole with mass Mbh can be cal- 
culated, if some parameters m, a, /3 and 6 are 
specified. The value of the vis cosity parameter 



NaravanI (|l996[ ) 



a is still a controversial issue, 
assumed a fairly large value, a = 1, in order to 
reproduce the maxim um luminosities o bserved in 
low-state XRBs, while lEsin et al.l (jl997l ) suggested 
a smaller value, a = 0.25, in their spectral model 
calculations. Recently, the influence of viscosity 
parameter a on the transition luminosity between 
the spectral states of accreting X-ray binaries was 
explored with a w ide range of a from 0.1 to 0.9 
doiaoTLi^llooi), which shows that the detected 
transition luminosity between the spectral states 
of many galactic accreting X-ray binaries can be 
fitted by tuning the value of a. In this work, we 
adopt two different values of a = 0.2 and 0.5 in 
the calculations, respectively. The parameter /3, 
defined as = -B^/Stt = (1 - l3)ptot, (ptot = 
Pg&s + Pm), describes the magnetic field strength 
of the accretion flow. We assume (3 — 0.8 in all the 
calculations. This parameter will mainly affect the 
radio spectrum from the source, while it has lit- 
tle effect on the X-ray emission, which we mostly 
focus on in this work. The parameter S describes 
the fraction of the viscously dissipated energy di- 



rectly going into electrons in the accretion flow. It 
was pointed out that a significant fraction of the 
viscously dissipated energy could go into electrons 
by magnetic reconn ection, if the magnetic fields in 
the flow are strong (jBisnovatvi-Kogan fc Lovelace 
I997LI2OOOI) . The value of 6 is still uncertain, and 
we adopt a conventional value of£ = 0.1 in all the 
calculations as that adopted in [Caol (|2007h . 

In this work, the global two-temperature struc- 
tures of the ADAF models are derived numeri- 
cally by solving the full relativistic hydrodynami- 
cal equati ons under t he mo deling techniques sug- 
gested bv iManmotol (|2000l ). The main improve- 
ments of this solution compar ed with th e previous 
works (e.g., G amm i c fc PophamI Il998 ) are that 
the advected fraction of the dissipated energy / in 
the energy equations is not treated as a parameter 
and the detailed radiation processes including syn- 
chrotron, comptonization, and bremsstr ahlung are 
calcul ated self-consistently. Although Manmotol 
(I2OOOI) bund that the solution with / = 1 is 
quite identical with the solution obtained by fully 
solving the energy equations, especially in the 
innermost part of the flow, there is a difference 
in the surface density in the outer regio n of the 
ADAF (see Figure 4 in iManmoto I2OOOI) . which 
is the dominant region contributing to the ther- 
mal X-ray iron line emissions we explor ed in this 
work. In Naravan fc Ravmond ( 1999f l. the ad- 
vected fraction of the dissipated energy / is as- 
sumed to be a constant i n solving the global 
structur es of the AD AFs fO.uataert fc Naraya: 
1999; G ammie fc Pop ham 1998; Esin ct al. 199 ' 
and the continuum emissions from the ADAFs 
are computed using the solved global structures 
of the ADAFs. Moreover, they have not calcu- 
lated the global structure of the ADAFs with 
winds, instead, they simply adopted the ra- 
dial velocity, angular velocity, and sound speed 
of the flo w derived from the global relativistic 
models of iGammie fc PophamI ( 19981 ) under the 
assumption of a constant mass accretion rate. 
The density of the accretion flow is then calcu- 
lated by applying the continuity equation with 
a variable r-dependent m ass accretion rate (see 
Quataert fc NaravanI 1999L for the details) . 



Given the values for parameters a, /3 and S, 
the mass accretion rate to can be tuned to fit 
the observed X-ray continuum spectrum when 
the outer radius of ADAF is fixed. In low-mass 
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BHXBs, where mass is transferred via Roche lobe 
overflow, the circularization radius of the incom- 
ing strea ni is large (flcirc ^ 10'' — 10^i?s, = 
2GM/c^) (|Esin et all Il997l ). We assume the di- 
mensionless outer radius of the ADAF surround- 
ing the black hole in GX 339-4 to be rout = 
-Rout/-Rs = 5 X 10'* in our calculations. A con- 
stant mass accretion rate m independent of ra- 
dius is adopted for the pure ADAF model. The 
jet/outflow models were widely used to account for 
the observed broadband spectra and the radio/X- 
ray correlations of GX 339-4 and other X-ray bi- 
naries in the low or off state, though the origin of 
the outfl ows in the ac c retion flows is still unclear 
Coron giu et al.l l2003t iMarkoff et aP l200ll 120031: 



It 



Tomsick et al. 2008f) . In this work, we employ a 



power-law r-dependent mass accretion rate, 



m 



t(r/ro 



(1) 



to describe the outflows from the ADAF, where 
the power-law index p represents the strength of 
the outflows in the accretion flow, and mout is the 
mass accretion rate at the outer radius. It can 
be easily found that the pure ADAF model is a 
speciflc case of ADAF with outflows model for p = 
0. The global structure of the ADAF is calculated, 
and then the spectrum is derived to reproduce the 
observed X-ray continuum by tuning the values of 
rhout and p. 

3. Calculations of line equivalent widths 

Based on the structure of the ADAF calculated 
as described in the last section, we can compute 
the thermal X-ray iron line emission from GX 339- 

4. The equivalent width of an emission line is de- 
fined as 

EWlinc = Liinc/L„{iyiinc), (2) 

where Lune is the total luminosity of the emission 
line, and ijy(i'une) is the spectral luminosity of the 
continuum at the energy of the line z^iinc- The total 
line luminosity inno is available by integrating the 
emission from the hot gas over the different shells 
of the accretion flow. 



-^line — 



nlir) eiir^ciTc) ^nrHdr Rl (3) 



where ric is the electron density at dimension- 
less radius r — R/Rg, r-m = Rm/Rs and rout = 



-Rout/^S are the inner and outer radii of the ac- 
cretion flow. The density Uc and the vertical half 
thickness H of the ADAF as functions of radius 
are given by the derived global ADAF structure. 
The quantity eunc is the line emissivity for unit 
density of plasma at a certain radius, which is a 
function of the local electron temperature. 

Given the electron temperature of the plasma, 
we use the standard software packag e Astrophys- 



ical P lasma Emission Code (APEC) (jSmith et al 



20011 ) to calculate the emissivity of the chosen 
line. We assume a solar abundance and ioniza- 



tion e quilibrium in the calculations. iPerna et al 



(|200rf ) :bund that departures from ionization equi- 
librium were small in a model with a^O.l, but 
they will be more severe in the models with larger 
values of a adopted here. The APEC code in- 
cludes coUisional excitation, recombination to ex- 
cited levels and dielectronic satellite lines. It ig- 
nores photoioniz ation, which is a few p ercent effect 
at most (Naravan fc Ravmond1ll999[ ). The spec- 
tral continuum luminosity Li,{viine) from the ac- 
cretion flow is calculated with the derived ADAF 
structure. We compute the He-like and H-like line 
equivalent widths of iron, EWi and EW2, centered 
at energy ~ 6.7 keV and ~ 6.97 keV, respectively. 
The solar metallicity Zq is adopted in the calcu- 
lations. The results can be scaled to other metal- 
licities Z using 



EWi,2(^) - ^EWi,2(Z0) 



(4) 



4. Results 



The system parameters (and companion star) 
of GX 339- 4 have not been directly measured. 
Hvnes et al. 1 (|2003h estimated a 1.7557 days or- 
bital period by analyzing the observed NIII 
(Bowen blend) lines and He II wings, and derived 
a minimum mass of 5.8Mq for the black hole in 
GX 339-4. The distance of GX 339-4 is estimated 
to be d > 6 kpc from the optical Na D line pro- 
files, and perhaps ev en to be as larg e as 1 5 kpc 



(|Hvnes et al.ll2004l) . IZdziarski et al.j (|2004l ) esti- 



mated the distance to be > 7 kpc based on Very 
Large Telescope (VLT) observation of GX 339-4. 
We adopt M = 10 Mq and d = 8 kpc in all the 
ca lculations, which are the same as those adopted 



Tomsick et al] (|2009l) 



As discussed in §2, most of the disk parameters 



4 



are fixed in the calculations of the global struc- 
ture of the ADAFs, except the viscosity parame- 
ter a, the parameter p describing the strength of 
the outflows, and the mass accretion rate at the 
outer radius of the ADAF rhout ■ The calculations 
are carried out with different viscosity parameters, 
a — 0.2 and a = 0.5, respectively. The calculated 
global structures of some ADAFs with different 
parameters are plotted in Figures 1 and 2. 

Based on the derived ADAF structur e, the con- 
tinuum spectrum can be calculated (see CaolbOOSl 
[2007, for the details). The value of the mass accre- 
tion rate at the outer radius of the ADAF mout and 
the outflow strength parameter p are tuned in or- 
der to fit the observ ed X-ray continuu m spectrum 
between 0.4-11 keV (|Gallo et al.ll2003l) . It is found 
that the values of rhout and p are degenerate, i.e., 
the X-ray continuum spectrum can be fairly well 
reproduced with different sets of the parameters 
rhout and p (see Figure [3]). 

Assuming a solar metallicity for the plasmas in 
the ADAF, we calculate the equivalent widths of 
the thermal He-like and H-like iron lines emitted 
from the accretion flows as described in §3. The 
results are listed in Table 1, where EWi(Z0) and 
EW2(^0) represent the equivalent widths of the 
thermal He-like and H-like iron lines, respectively. 
The adopted parameters are also shown in Ta- 
ble 1. The total equivalent widths of the He-like 
and H-like iron lines emitted from the ADAFs sur- 
rounding the black hole in GX 339-4 are plotted 
in Figure 4 as a function of the outer mass accre- 
tion rate for the models with different strength of 
the outflow. The results show that the predicted 
EWs of the He-like and H-like iron lines from the 
pure ADAFs are too small to be observed. For 
the ADAFs with outflows, the EWs of the thermal 
iron lines significantly increase with the strength 
of the outflow. 

5. Discussion 

In this paper we have explored the thermal X- 
ray line emission from the accretion flow surround- 
ing the black hole in the galactic binary GX 339-4. 
We construct the ADAF models for the observed 
X-ray spectrum of GX 339-4 with different param- 
eters including the viscosity parameter a, the mass 
accretion rate rh and the outflow strength p. We 
compute the equivalent widths of the helium-like 




R/R, 



Fig. 1. — Electron density distributions of the 
ADAF models with a = 0.2 and a = 0.5, re- 
spectively. The black solid line corresponds to the 
pure ADAF model with parameters a = 0.2, m = 
0.001. The red and blue solid lines correspond to 
the ADAF with outflows models with parameters 
a = 0.2, rhout = 0.008, p = 0.25, and a = 0.2, 
rhout = 0.025, p — 0.45, respectively. The black 
dot-dashed line corresponds to the pure ADAF 
model with parameters a = 0.5, rh = 0.001. The 
red and blue dot-dashed lines correspond to the 
ADAF with outflows models assuming parameters 
a = 0.5, rhout = 0.01, p = 0.25, and a = 0.5, 
rhout = 0.045, p = 0.45, respectively. 




10 10 
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Fig. 2. — Same as Figure 1, but for the electron 
temperature distributions. 
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Table 1 

Equivalent widths EWi(Z0) and EW2(^0) in units of eV of the thermal helium-like and 
hydrogen-like x-ray iron lines from the adafs 



a 


TOout 


P 


EWi(cV) 


EW2(cV) 


0.2 


0.001 


0. 


0.7 


0.7 


0.2 


0.008 


0.25 


28.9 


24.3 


0.2 


0.01 


0.3 


54.5 


41.1 


0.2 


0.015 


0.35 


103.9 


69.0 


0.2 


0.02 


0.4 


196.2 


110.1 


0.2 


0.025 


0.45 


368.0 


169.2 


0.5 


0.001 


0. 


0.4 


0.5 


0.5 


0.01 


0.25 


7.4 


6.5 


0.5 


0.015 


0.3 


17.0 


13.6 


0.5 


0.025 


0.35 


38.8 


27.5 


0.5 


0.035 


0.4 


77.2 


48.2 


0.5 


0.045 


0.45 


141.0 


77.1 



_ 10'"- 
'« ^ ■ 

> 

_l 

10"- 




PjO.45 

p=0.40 
O 

p=D.35 



p=0.30 
O 

p=0.25 



10' 



10" 



10' 



v(Hz) 



Fig. 3. — Spectra of ADAF models with a = 0.2 
and a = 0.5, respectively. The black dotted line 
corresponds to the pure ADAF model with param- 
eters a = 0.2, TO — 0.001. The red and blue dot- 
ted lines correspond to the ADAF with outflows 
models with parameters a = 0.2, rliout = 0.008, 
p = 0.25, and a = 0.2, TOout = 0.025, p = 0.45, re- 
spectively. The black dot-dashed line corresponds 
to the pure ADAF model with parameters a = 0.5, 
m = 0.001. The red and blue dot-dashed lines 
correspond to the ADAF with outflows models as- 
suming parameters a = 0.5, mout — 0.01, p = 
0.25, and a = 0.5, TOout — 0.045, p = 0.45, respec- 
tively. The black tie repres ents the obs e rved .4-11 
keV X-ray spectrum from i Gallo et al.l (|2003l) . 



o 
+ 



10"' 



Fig. 4. — Total equivalent widths of the He-like 
and H-like iron lines from the ADAFs as a func- 
tion of the outer mass accretion rate. Circles and 
stars represent the data for models with a = 0.2 
and a = 0.5, respectively. The power-law indexes, 
p, in Equation ([l} are indicated in the figure for 
different ADAF models, while p = {) corresponds 
to the pure ADAF models. 
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and hydrogen-like thermal iron lines from the ac- 
cretion flows. 

From Table 1 and Figure 4, we can see that 
the thermal X-ray line emission from the ADAFs 
without outflows is always much lower than emis- 
sion from the ADAFs with outflows. Moreover, 
the equivalent widths of the He-like and H-like 
thermal iron lines significantly increase with the 
mass accretion rate rriout at the outer radius of 
the ADAF and the outflow strength parameter 
p. In fact, these two parameters are not inde- 
pendent, i.e., the outflow strength parameter p is 
required to increase with the mass accretion rate 
rhout at the outer radius to reproduce the observed 
X-ray continuum spectrum. The sensitive depen- 
dence of EWi and EW2 on the two parameters 
can be explained by the very different density pro- 
files (see Figure 1), because the temperature pro- 
files of the plasmas are quite similar for ADAFs 
with different parameter s adopted (see Figure 2, 
and iNaravan fc Yil 119951) . As showed in Equa- 
tion ([2]), the equivalent width of a thermal line 
is the ratio of the line luminosity and the spec- 
tral luminosity at the energy of the line. All our 
derived spectra are required to fit the observed 
X-ray continuum at ~ 6.83 keV. Thus, the equiv- 
alent widths of thermal iron lines will mainly de- 
pend on the line luminosities. From Equation ([3l). 
the line luminosity emitted from a certain annulus 
of the accretion flow is determined by the density 
and the line emissivity, which is a function of the 
electron temperature. As the temperature pro- 
files are quite similar for the models with different 
parameters, the line luminosity is mostly deter- 
mined by the density distribution of the gas in 
the accretion flow. For the models with larger p, 
i.e., stronger outflows, the density profiles are fiat- 
ter, which means that much more gases are blown 
away from the accretion fiow. This leads to higher 
electron densities in the outer region of the accre- 
tion flow (r > 10'^) with the electron temperature 
< 10* K, where the thermal He-like and H-like 
iron lines emitted. Therefore, the line luminosity 
and the equivalent widths of the thermal iron lines 
are higher in these models. 

The change of the equivalent widths of the ther- 
mal iron lines emitted from the ADAFs with the 
value of the viscosity parameter a is explored in 
this work. Comparing the two groups of the line 
equivalent widths for different values of viscosity. 



a — 0.2 and a — 0.5, we flnd that the EWs of 
the lines from the ADAFs with a = 0.2 are always 
larger than those from the ADAFs with a = 0.5 if 
the same value of p, i.e., strength of the outflows, 
is adopted. This is because the larger the viscos- 
ity is, the quicker the gas flows in, which leads to 
smaller densities of the gas in the outer region of 
the ADAF (see Figure 1). 

We theoretically conclude in this work that the 
galactic X-ray binary GX 339-4 may produce X- 
ray iron lines with very large equivalent width in 
the quiescence. For reasonable choice of model pa- 
rameters, e.g., a = 0.2, TTiout = 0.25 and p — 0.45, 
the total equivalent width of the He-like and H-like 
iron lines from the ADAFs with moderately strong 
outflows reaches to > 500 eV, when the 0.4-11 
keV X-ray flux is about 3.8 x 10~^'^erg cm~^ s~^. 
If these strong thermal lines are really observed, 
it may give a clue to the accretion mode of the 
source in the off state and provide strong evi- 
dence t hat the outflows/w inds exist in GX 339-4. 



In fact, iFeng et al.l ()2001l) showed that a promi- 
nent thermal X-ray iron line, with the equivalent 
width of 595^^72 and the line energy centered 
at 6.83tj^:i^ keV, is present in the RXTE (PGA) 
observation of GX 339-4 in 1999 when the source is 
in the off state with the integrated 3-20 keV flux 
of 7 X lO^^^crg cm^^ s^^. However, we need to 
bear in mind that the observation with RXTE may 
include possible contamination from the Galactic 
ridge or other sources. Observations with ASGA, 
Suzaku and Ghandra have detected an iron line 
emission at 6.7 keV and suggested it to be emit- 
ted by the diffuse hot plas ma and/or produced 
by p oint sources (se e, e.g., Kovama e t al. 1986t 

2.009). T hus, we 



py 1 
Iran 



Tanak a ,2002; Revniv tsev et al 



select the Ghandra observation (jGallo et al.ll2003[ ) 
as the quiescence flux in our calculations. 

Due to the different sources we explored and 
different parameters we used in the calculations, 
it is not easy to directly compare the equivalent 
widths of thermal iron lines from GX 339-4 derived 
in this work with those fr o m V4 04 Gyg predicted 
by iNaravan fc Ravmond ( 19991 ). Therefore, we 
recalculate the equivalent widths of the thermal 
iron lines emitted from V404 Gyg in quiescence 
with our ADAF mo del for the same model parame - 



ters as those used in lNaravan fc Ravmond l (ll999l) 



The comparison of the resu lts is shown in Table 2 
The value of parameter /? in lNaravan fc Ravmond 
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( 19991 ) (defined as /3 = Pgas/pmag) has been trans- 



lated to the value with the definition in this work. 
T he equivalent width of the He- like iron line EWi 
of Naravan fc Raymond ( 1999t ) is the sum of the 
several He-like lines with the energy between 6.65 
keV and 6.75 KeV which can be compared with 
the result of this work, because we have already 
combined the equivalent widths of He-like lines 
near the energy of 6.7 keV in bins of 100 eV 
width which we believe to be more comparable 
with the observation. Moreover, we note that 
the equivalent width o f the H-like iron line of 
Naravan fc Ravmond ( 1999^ only includes one 
line Fe XXVI A1.780 while the data of this work 
is also the sum of the H-like lines in bins of 100 
eV width, i.e. with the energy between 6.92 
keV and 7.02 keV. The emission in the 100 eV 
bin is larger than tha t of th e line calculated by 



Naravan &: Ravmond 



(|l999l ). The calculated X- 
ray continuum derived from the approach of this 
work is a bit l arger than the ASCA observation 



of V404 Cyg (iNaravan et al.l Il997l) if the same 



mass a ccretion rate mout as lNaravan fc Ravmond 
(|1999[ ) is adopted. Using the same model pa- 
rame ters 6,a,0 and the s treng th of the wind p 



as iNaravan fc Ravmond I ()1999D . our calculation 
shows that the best fitted mass accretion rate is 
Wout — 0.004 (the results are listed in Table 2). 
The discrepancy of the resu l ts bet ween this work 
and Naravan fc Ravmond ( 19991 ) is significant, 
which is mainly due to the difference of ADAF 
structure calculations. 

We thank the anonymous referee for very 
helpful suggestions. This work is supported 
by the NSFC (grants 10778621, 10703003, and 
11078014). 
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this work 
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0.006 
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this work 
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0.4 


0.004 


122 


79 
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